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We report on a study of heavy quark bound states containing an additional excitation of the 
gluonic degrees of freedom. To this end we employ the NRQCD approach on coarse and asymmetric 
lattices, where we discard vacuum polarisation effects and neglect all spin-correction terms. We find 
a clear hybrid signal on all our lattices (a^ = 0.15 . . . 0.47 fm). We have studied in detail the lattice 
spacing artefacts, finite volume effects and mass dependence. Within the above approximations we 
predict the lowest lying hybrid excitation in Charmonium to be 1.323(13) GeV above the ground 
■ state, where we use the IP-IS splitting to set the scale. The bottomonium hybrid was found to be 

^ ' 1.542(8) GeV above its ground state. 
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Gluonic excitations are ideal objects for investigating the nonperturbative nature of the gluonic degrees of freedom 
in QCD. Hybrid mesons can be thought of as hadronic bound states with an additional excitation of the gluon flux. 
This can give rise to states with non-conventional quantum numbers and has triggered an intense experimental and 
• theoretical search for such particles. Previous predictions for the energies of hybrid states come from phenomenological 
CN models static potential models and lattice simulations with propagating quarks [||-||. So far, lattice QCD is 
the only approach in which hybrids can be treated from first principles. However, the errors from such calculations 
on isotropic lattices are still much larger than for conventional states. This is because the correlation functions decay 
too rapidly when the excitation energy is very large with respect to the inverse lattice spacing. To obtain a similar 
\ signal-to-noise ratio one needs a much finer resolution in the temporal direction. 

While in recent years there has been much progress in obtaining more reliable results from improved actions on 
CIhI spatially coarse lattices, it has also been demonstrated that anisotropic lattices can be employed to accommodate 
^ • different physical scales on the same lattice. In particular the study of glueball states on coarse and anisotropic lattices 
has prompted us to study heavy hybrid states on such lattices in order to increase both the scope and the precision 
of a previous calculation Q significantly. There a non-relativistic approach (NRQCD) was used for the heavy b quarks 
on an isotropic lattice with a « 0.08 fm. NRQCD has frequently been employed to allow high precision measurements 
for the 66-system . Also the combination of improved gluon actions and the NRQCD approach for heavy quarks 
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I has already been used to determine the spectrum of conventional quarkonia |13-15f|. Previous attempts to measure 

heavy hybrid states on the lattice were reviewed in p^ . An initial preliminary study using asymmetric lattices for 
Bottomonium hybrids was presented in P,[l7|. 

In this paper we implemented such an efficient approach to study in detail the lattice artefacts and finite size effects 
for both the bbg and ccg hybrid. For the latter state we obtain excellent agreement with a relativistic simulation on 
isotropic lattices [||. As with this other study we have neglected dynamical sea quark effects, but we have succeeded 
in lowering the statistical error to about 1%. 

In our study we generated the gauge field configurations using a tadpole-improved action which has been employed 
by different groups |@,|l5|: 

a:,i>j x,i 

Here Pjj and i?ij denote the trace of the standard spatial plaquette and rectangle, respectively. Where the index t 
appears the plaquette/rectangle extends only one link into the temporal direction. This theory has two parameters, 
P and ^, the second of which determines the asymmetry of our lattices. At tree level the "aspect ratio" is ^ = ag/at, 
where and at are the spatial and temporal lattice spacings, respectively. From 0| we note that the radiative 
corrections to this relation are small when tadpole improvement is implemented, as described below. 
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The action in Equation ^ is designed to be accurate up to ©(a^, aj) classically. To account for radiative corrections 
all spatial gauge links, Ui, are self-consistently tadpole-improved with = (0|(l/3) Pij|0)^/'', as suggested in [Q. Such 
a mean-field treatment was demonstrated to reduce significantly the leading corrections, which are due to unphysical 
tadpoles in lattice perturbation theory. Our scaling analysis shows that errors 0{aa^) are indeed negligible if the 
lattice spacing is sufficiently small. Since we will use this action only for small temporal lattice spacings we expect 
Oi^aaj) errors to be very small. Therefore we have not employed the tadpole improvement description for the temporal 
gauge links. Our results for two different aspect ratios ^ justify such an assumption. 

To propagate the heavy quarks through the lattice we expand the NRQCD Hamiltonian correct to 0{mv^), where 
all spin-dependent terms are absent. This accuracy was already employed in The only additional improvement is 
to correct for temporal and spatial lattice spacing errors by adding two extra terms (cy, cg) to the evolution equation 
in0]: 

In those operators all spatial links are also tadpole-improved using the same Us as for generating the configurations. 
After this modification we take the tree-level values for all the coefficients in the Hamiltonian. In this case, C7 = cg = 1. 

For the non-relativistic meson operators we only used the simplest possible choices and have not tuned the overlap 
to optimise the signal. The gauge-invariant construction of S-state and P-state operators is described in For 
the magnetic hybrid signal studied here, we have inserted the lattice version of the magnetic field into the QQ'State 

^H,ix)^^P\x)B,x{x) , ''H,k{x)^^P\x)ajBkxix) . (3) 

For the leading order in the NRQCD Hamiltonian, the operators in Equation || create a whole set of degenerate states: 
1 , ^ , 1 ^ ^ 2 ^ . These are the spin-singlet and spin-triplet states with zero orbital angular momentum, including 
the exotic combination 1 ^ . This degeneracy will be lifted when higher order relativistic corrections are re-introduced 
into the NRQCD Hamiltonian. We expect this to be a small effect, in as much as the heavy quarks are very slow 
in the shallow hybrid potential By the same argument, we expect hybrid states with additional orbital angular 
momentum to be almost degenerate as it was observed in . The definition of the operators in Equation ^ has been 
augmented by a combination of fuzzing p9| for the links and Jacobi-smearing for the quark fields [ pO[ . No effort has 
been made to optimise the signal further, but one could do so if even higher precision is needed, or if higher excited 
states are to be determined. To extract the hadron masses we simply fit the meson correlators to a single exponential; 

a(i) = (i/IWi?a(0))=A„e-""* . (4) 

We measured correlators every 10 trajectories in the Monte Carlo update, and for the error analysis we binned 50 
such measurements into one. After this binning we still have an ensemble of 100-1000 configurations depending on 
the lattice and the state of interest. As it can be seen from the representative example of an effective massplots in 
Figure 0, the data is very good and the goodness of the single exponential fits is always bigger than Q = 0.1, which 
we called acceptable. In Tables ^ and R we present our results and the simulation parameters. 

To determine the lattice spacing , we used the IP-IS splitting in Charmonium and Bottomonium. As expected, 
the values from Charmonium are smaller than those from Bottomonium, because in the quenched approximation the 
coupling does not run as in full QCD. At (/3,^) = (2.4,5) we observe a 16% effect. In order to give another estimate 
of quenching errors we also determined the radial excitations, nS^ and calculated the ratio Rsp= 2S-1S/1P-1S, which 
can be compared with the experimental value of 1.28 for Bottomonium. At (/3,C) = (2.7,5) we find i?sp — 1.424(89). 
From these findings we quote quenching errors of 10 — 20%. Several suggestions have been made how to measure the 
spatial lattice spacing [p^ , ^ , but to convert our results into dimensionful numbers we can use a^^ throughout. 

We have also tested the velocity expansion and included relativistic corrections up to 0{mv^) for some of our 
lattices. At this level of accuracy we could not resolve any significant change in our results. A more detailed analysis 
of the spin structure is subject of a future project. 

In this study we were mainly interested in the gluonic excitations of heavy quarkonia. For this purpose it was 
irrelevant to adjust the quark masses to their exact values. For some lattices we have changed the quark masses by 
25% and did not find any noticeable change in the ratio i?H = (1H-1S)/(1P-1S). 

Finite volume effects were a source of immediate concern for us. This is because hybrid states are expected to 
reside in a very fiat potential |^ . The bag model also suggests a very large bound state as the result of the gluonic 
excitation [^. As shown in Figure |2| we have found that for spatial extents of 1.2 fm or larger the masses of all 
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Bottomonium and Charmonium states remain constant within small statistical errors. For even smaller volumes we 
can resolve a slight increase in the mass of the bbg hybrid. In other words, the spatial extent of the hybrid excitation 
seems to be almost independent of the heavy quark mass. In this sense the bbg hybrid is more difhcult to calculate, 
as we need similar volumes but finer lattices than for Charmonium. 

Finally, we carried out a scaling analysis to demonstrate that discretisation errors are under control. This is of 
utmost importance for the NRQCD approach since one cannot extrapolate to zero lattice spacing in this effective field 
theory and one has to model continuum behaviour already for finite lattice spacings. In Figure ^ we show the scaling 
of the hybrid excitation above the ground state. From this, one can see that we have found convincing scaling windows 
for both Bottomonium and Charmonium. Scaling violations can only be seen on the coarsest lattices (a^ > 0.36 fm 
for ccg and Ug > 0.19 fm for bbg), but this is not totally unexpected - it is questionable how well our simple minded 
implementation of the tadpole prescription works to remove the 0{ap'^a'^) errors for heavy quark systems on such 
coarse lattices. In the case of bbg we also plot the results for two different aspect ratios. Both results are consistent 
and confirm the initial assumption of small temporal lattice spacing errors. To quote our final result for the lowest 
lying hybrid excitations we take the averaged value of all the results in the scaling region and employ the experimental 
values for the IP-IS splitting to set the scale. We find 1.323(13) GeV for the case of Charmonium and 1.542(8) GeV 
for the first gluonic excitation in Bottomonium, in good agreement with a previous estimate of 1.68(10) GeV 

In conclusion, we have demonstrated the usefulness of coarse and anisotropic lattices for the nonperturbative study 
of gluonic excitations in heavy quark systems. Furthermore, this should also be considered a success of the NRQCD 
approach, which allowed us to predict the lowest lying Charmonium hybrid state at the same mass as from a relativistic 
calculation [Q. Apart from the very accurate predictions for hybrid quarkonia, it is also interesting to notice that 
all of the above results could be obtained in a comparatively short period of time. Whereas our present calculation 
confirms that the bbg hybrid will lie above the S+S threshold for decay into B-mesons, the issue whether it may be 
found below the S+P threshold has to be decided in a simulation where dynamical sea quarks are included in order 
to control this last remaining systematic error. 
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at the University of Tsukuba. This work is supported in part by the Grants-in-Aid of Ministry of Education (No. 
09304029). TM, HPS and AAK are supported by the JSPS Research for Future Program, and SE and KN are JSPS 
Research Fellows. 



[1] N. Isgur,J. Paton, Phys.Rev. D31, (1985), 2910. P. Hasenfratz et al, Phys.Lett. B95, (1981), 299. T. Barnes et al, 
Phys.Rev. D52, (1995), 5242. 

L.A. Griffiths et al, Phys.Lett. B129, (1983), 351. S. Perantonis, C. Michael, Nucl. Phys. B347, (1990), 854. 
K.J. Juge et al, Nucl.Phys.(Proc.Suppl.) 63 (1998), 326. 
C. Bernard et. al, Phys.Rev.D56 (1997), 7039. 
P. Lacock et. al, Phys.Rev.D54 (1996), 6997. 

T. Manke et al, Nucl.Phys.(Proc. Suppl.) 63 1998, 332; S. Collins et al, Nucl.Phys.(Proc. Suppl.) 63 1998, 335. 
C. Morningstar, M. Peardon, Phys.Rev.D56 (1997), 4043. 
T. Manke et al, Phys.Rev. D57 (1998), 3829. 
C.T.H. Davies et al Phys.Rev. D50 (1994), 6963. 
T. Manke et al, Phys.Lett. B408 (1997), 308. 
N. Eicker et al, Phys.Rev. D57 (1998), 4080. 
C.T.H. Davies et al, Phys.Lett. B345 (1995), 42. 
H.D. Trottier, Phys.Rev. D55 (1997), 6844. 
M. Alford et al, Phys. Lett. B361 (1995 ), 87. 

LT. Dru mmond et al, ^ep-lat/980917C| , Nucl.Phys. (Proc. Suppl.) m press. 



[iep-lat/981102l|, Nucl.Ph ys.(Proc.Suppl.) in press 



J. Kuti, 

K.J. Juge et a/., |hep-lat/9809098|, Nucl.Phys. (Proc. Suppl.) in press. 



CP. Lepage, P.B. Mackenzie: Phys.Rev. D48 (1993), 2250. 

M. Albanese et al, Phys.Lett. B192 (1987), 163. 

S. Giisken et al, Nucl.Phys. (Proc.Suppl.) (1990), 361. 

T.R. Klassen, Nucl.Phys.B533 (1998), 557. 

K.J. Juge et al, Nucl.Phys.(Proc.Suppl.) (1998), 543. 



3 



(0 £) 
Volume 


(1 7 5) 
4^ X 40 


(19 5) 
4^ X 40 


(2 2 5) 
8^ X 40 


(2 4 5) 
8^ X 40 


(2 4 5) 
8^ X 40 


Us 


0.7370 


0.7568 


0.7841 


0.7997 


0.7997 


asTTlc 


3.0 


2.66 


2.0 


1.62 


1.62, 0{mv'^) 


P-S 

a^VGoV 


0.2196(18) 
2.084(18) 


0.1689(26) 
2.709(42) 


0.1299(13) 
3.522(37) 


0.1068(21) 
4.286(86) 


0.1047(42) 
4.37(18) 


H-S 

H-S/P-S 
H-S/GcV 


0.5713(98) 

2.602(49) 

1.191(23) 


0.4860(48) 

2.877(53) 

1.317(24) 


0.3821(33) 

2.928(53) 

1.346(18) 


0.3048(17) 

2.855(59) 

1.306(27) 


0.3011(24) 

2.87(12) 

1.316(54) 



TABLE I. Results for Charmonium. The dimensionful numbers in the scaling region are given in boldfaced characters. From 
their average we obtain 1.323(13) GeV for the lowest lying hybrid excitation from our sinmlation with accuracy 0{mv^, a^, of). 
In the last column we give the spin averaged results from a higher order accuracy 0{mv^,at,at)- 



Volume 

Us 

asTUb 


(2.4,5) 
8^ X 40 
0.7997 
4.73 


(2.6,5) 
6^ X 40 
0.8139 
3.50 


(2.7,5) 
9^ X 40 
0.8193 
3.15 


(2.7,5) 

9'^ X 40 
0.8193 

3.15, 0{niv'') 


(2.5,3) 
6^ X 40 
0.8100 
4.10 


(2.6,3) 
8^ X 40 
0.8193 
4.00 


(2.8,3) 
lO-'* X 40 
0.8314 
3.33 


P- S 

a^VGeV 


0.08665(68) 
5.075(43) 


0.07444(65) 
5.907(54) 


0.06573(62) 
6.690(67) 


0.0636(23) 
6.91(25) 


0.1388(10) 
3.169(25) 


0.13068(82) 
3.365(23) 


0.10988(77) 
4.002(30) 


H-S 

H-S/P-S 
H-S/GeV 


0.30988(84) 

3.576(30) 

1.573(14) 


0.2591(11) 

3.481(34) 

1.531(15) 


0.2275(36) 

3.462(65) 

1.522(29) 


0.2305(46) 

3.62(15) 

1.594(67) 


0.4882(41) 

3.517(39) 

1.547(18) 


0.4568(21) 

3.496(27) 

1.537(13) 


0.3852(32) 

3.506(38) 

1.542(17) 



TABLE H. Results for Bottomonium. Prom the average over the scaling region we obtain 1.542(8) GeV for the lowest lying 
bbg hybrid, when the P-S splitting is used to set the scale. In column 5 we show results with accuracy 0{mv^,at,at). 
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FIG. 1. A representative effective massplot for the S, P and 
hybrid state in Bottomonium at (/?, ^) = (2.7, 5), asTUb = 3.15 
on a 9^ X 40 lattice. 
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FIG. 2. Finite volume analysis for Charmonium and Bot- 
tomonium. Wc plot the dimonsionloss energies for different 
states against the inverse spatial extent, 1/L, of the lattice. 
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FIG. 3. Scaling of the hybrid excitation, H — S. We plot 
the ratio Rn = {H — S)/{P — S) against the squared spatial 
lattice spacing; Os = ^at at tree level. We also show the 
results from a previous calculation [8] on a symmetric lattice 
at P — 6.0 (burst). To display the result for the lowest ccg 
hybrid at /? = 6.15 [4], we use their value 1.32(8) GeV and 
the experimental P-S splitting in Charmonium (triangle). 
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